Theoretical noise analysis on a position-sensitive Metallic Magnetic 

Calorimeter 


Stephen J. Smith’ 

NASA Postdoctoral Program Fellow, NASA Goddard Space Flight Center, Code 662, Greenbelt, MD 20771, USA 

Abstract 

We report on the theoretical noise analysis for a position-sensitive Metallic Magnetic Calorimter (MMC), 
consisting of MMC read-out at both ends of a large X-ray absorber. Such devices are under consideration as 
alternatives to other cryogenic technologies for future X-ray astronomy missions. We use a finite-element 
model (FEM) to numerically calculate the signal and noise response at the detector outputs and investigate 
the correlations between the noise measured at each MMC coupled by the absorber. We then calculate, using 
the optimal filter concept, the theoretical energy and position resolution across the detector and discuss the 
trade-offs involved in optimising the detector design for energy resolution, position resolution and count rate. 
The results show, theoretically, the position-sensitive MMC concept offers impressive spectral and spatial 
resolving capabilities compared to pixel arrays and similar position-sensitive cryogenic technologies using 
Transition Edge Sensor (TES) read-out. 
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1.0 Introduction 

Cryogenically cooled detectors such as Transition Edge Sensors [1], Metallic Magnetic Calorimeters [2] and 
Superconducting Tunnel Junctions [3], are under development by numerous groups worldwide for 
applications in X-ray astronomy. Large format pixel arrays (typically 32x32) [4] are required for future 
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astronomy missions such as NASA’s Constellation-X [5], The goal is to achieve full width at half maximum 
(FWHM) energy resolution of a few eV in the soft X-ray energy range. 

Large absorber devices, using microcalorimeter or STJ read-out at both ends of a linear or segmented 
absorber [6-8], offer alternatives in achieving wide field coverage compared to pixel arrays. In the TES 
Distributed Read-Out Imaging Device (DROID), the sum of the two TES signals gives the energy of the 
absorbed photon and the difference between the signals gives the position information. Further details on 
these devices are reported elsewhere [8-10]. The development of such detectors is driven by the desire to 
achieve wide field coverage with a reduced number of read-out channels compared to the pixel arrays, thus 
providing a technologically simpler alternative and, the capability to extend the focal plane area for missions 
such as Constellation-X. 

The position-sensing concept used in the TES DROID can, in principle, be implemented using MMC 
technology. In this paper we examine the theoretical performance of the MMC DROID concept. 

An MMC consists of an X-ray absorber thermally coupled to a material whose magnetisation obeys a 1/T 
Curie law dependence, at temperatures of the order 50 mK. The absorption of energy 5E in a detector with 
volume V and heat capacity C causes a change in magnetic moment 8m, proportional to the deposited energy. 
The sensor is magnetised by a small B-field typically of the order a few mT. Changes in magnetic flux 
caused by X-ray absorption are measured using a SQUID (Superconducting Quantum Interference Device) 
magnetometer, which is directly flux coupled, or transformer coupled, to the sample. The magnetisation in 
MMCs comes from a dilute concentration of paramagnetic ions. This dilute concentration of randomly 
distributed ions in a host material is essential to limit the spin-spin interaction between ions, which result in a 
decrease in sensitivity. Unlike a dissipative thermistor such as a TES, MMCs require no bias power therefore 
dissipate no heat into the detector substrate and require no active bias circuitry. Thus the MMC is not subject 
to the same electrical and thermal design constraints as a thermistor and potentially, large format array 
development is technologically less challenging. Early results in single pixel MMCs have already 



demonstrated comparable resolution to other calorimeter types, A E = 3.4 eV FWHM [2] for 5.9 keV X-rays. 
Yet further improvements in resolution are theoretically predicted for fully optimised devices. MMCs 
therefore offer an exciting prospect for focal plane X-ray detectors. 

2.0 Detector signal and noise response 

We have established a linear finite-element model (FEM) to estimate the signal and noise response of the 
MMC DROID concept, similar to those described in refs. [1 1,12] for position-sensitive TES detectors. In the 
FEM, the MMCs are represented as each having two subsystems. Firstly, the electronic system of the MMC, 
including the embedded paramagnetic ions and the host absorber material, with a heat capacity C e and 
temperature T e . Secondly, the MMC spin system with heat capacity C z and temperature T z . The thermal 
conductance between the spin system and the electron system is G z and the thermal conductance between the 
electronic system and the bath, at temperature 7*, is G. The electron system is then connected to an absorber 
of length L, consisting of a series of N (typically N > 100) absorber heat capacities with overall thermal 
conductance G a bs and heat capacity C a b S - Figure 1 shows a schematic diagram of the FEM representing the 
MMC detector. 

The noise sources in the detector come about from thermodynamic fluctuations between the various sub- 
systems. These are: thermal fluctuation noise between the MMC and bath P„, thermal fluctuation noise in 
each conductive element between absorber nodes P^tbs and thermal fluctuation noise between the MMC spin 
system and the electron system P z . All these sources of noise have an intrinsic spectral density in the form 

yj4k h T 2 G [13], where kb is Boltzmann’s constant, G is the thermal conductance between the relevant sub- 
systems and T is the temperature. The time evolution of the detector can then be represented by a series of 
coupled, linear, differential equations. Following the methodology of ref. [12] for a similar TES based 
DROID, this system of differential equations can be generalised as a single matrix equation and solved to 
find the signal and noise response at the detector outputs. The total noise, in terms of temperature 


fluctuations in the spin system, when the MMC signals are summed, can be written (in the frequency 
domain) as: 


(K 2 (/)f) = (hoof) + (h 2 (/)| 2 ) + 2R e (i\^*(/)) , i 

where (|./Vi(/)| 2 ) and ( jA 7 " 2 (/)| 2 ) are the noise powers in each MMC spin system and (|A r ijV 2 *(/)|) is the noise 
covariance, which arises due to the coupling of the sensors via the absorber. 

In this analysis we consider a simple geometry consisting of a MMC detector with a circular pickup coil-up, 
which is then transformer coupled to a SQUID. We note that different flux coupling geometries are also 
possible and yield slightly different optimisations conditions [14]. 

Considering a simple system of non-interacting spins, the flux change <50 in a circular pickup loop of radius 
r, due to a change in temperature 8T is given by: 
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where Gk is a dimensionless factor which depends on the geometry between the spins and the pickup loop. 
The change in magnetisation with temperature dM/dT, is given by the thermodynamic expression [15]: 


dM C z 1 
dT ~ V B’ 


where V is the sensor volume and B is the applied magnetic field. In reality the situation is more complex and 
an exchange interaction will exist between the spins. The existence of - Ruderman Kittel Kasuya Yosida 
(RKKY) [16] - exchange interactions will have a detrimental effect on MMC performance by reducing 


dM/dT and hence, the device sensitivity [15]. In this regime: dM/dT & x g 2 B T 2 , where x is the paramagnetic 
ion concentration, g is the g-factor of the spins, B is the magnetic field and T is the temperature. A large 
dM/dT is essential to ensure that the performance of the device is not limited by the flux noise from the 
SQUID. This requires careful optimisation of the device design and operating parameters, such as x and B, to 
optimise the MMC resolution for a particular flux coupling scheme. Assuming a transformer coupled MMC, 
the value of dM/dT to optimise the energy sensitivity, is numerically shown to be proportional to g' l a\ 
where a is the magnitude to RKKY exchange interactions (for a typical Er:Au system g = 6.8 and a = 5) 
[16]. 


Using Equation 2 the intrinsic detector noise can be converted into a flux noise in the pick-up coil: 


(I^mmcI 2 ) = (^j f(K(/)| 2 ) + (|^ 2 (/)| 2 ) + 2Re(jV, N* 2 (f)) 
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In addition to the intrinsic detector noise we must also consider noise from the SQUID read-out chain. We 
assume the pickup coil is flux coupled to the SQUID, therefore changes in flux in the pickup coil 5<& ex , 
generate a change in flux in the SQUID <50 s given by: 


<50 = ~ - /;t - <50 
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where L p is the inductance of the pickup coil, L, is the inductance of the input coil and k is a coupling 
constant close to unity, M ls is the mutual inductance between the SQUID and the input coil. Thus the flux 
noise referred to the pickup coil is: 
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Using typical SQUID parameters (L, = L p = 50 pH and L s = 30 pH), and assuming an intrinsic SQUID noise 
of O s = 0.2 pdVVlHz, we calculate a flux noise in the pickup coil of O m = 0.5 pOo/VHz. 

Magnetic Johnson noise will also be a contributing factor [15]. Magnetic Johnson noise is a flux noise 
caused by the thermal motion of charge carriers in the sensor, the absorber and other metallic materials in 
close proximity to the sensor. Using the data from ref. [14] for a typical Er doped Au sensor, the magnetic 
Johnson noise will be ~ 0.1 pOo/VHz. The metallic absorber is also a potential source of magnetic 

Johnson noise. Again, from ref. [13] the flux noise <f>j a bs is estimated to be ~ 0.2 p<J> 0 A/Hz. 

The noise from the two SQUIDs and sources of magnetic Johnson noise will add incoherently with the 
detector noise, therefore, the total noise is: 


M = (Kmc(/)| 2 ) + 2(|<M 2 ) + 2(|<D jmmc | 2 ) + 2 ^ Jabt \ 2 ) . 7 

Figure 2 shows, using the device parameters in Table 1, the calculated flux noise spectral density in the 
pickup loop for one of the MMCs, <£„i. Figure 3 then shows the total noise spectral density $>mmc when the 
signals from both MMCs are summed together. To examine the intrinsic detector limits we set the magnetic 
Johnson noise and the SQUID noise initially to zero. 

Modelling of the noise spectrum reveals the frequency dependent correlations between the individual noise 
sources when sampled on each MMC. This is an identical situation to TES DROID except we have different 
sources of detector noise [10]. Both MMCs will sense all noise contributions via the coupling of the 


absorber; however, one MMC will sense noise attenuated and phase shifted relative to the other. When the 
two signals are summed, a frequency dependent correlation will exist. This is most clearly indicated in the 
absorber phonon noise, where at low frequency, the noise contributions tend to cancel. These correlations are 
quantified from the covariance (<I>i<I>2*(/)) of the data (see Figure 4). 

The phonon noise in the absorber is heat moving between adjacent heat capacities. At low frequency one 
MMC will detect a drop in temperature while the other detects an increase, thus when the signals are added, 
the noise effectively cancels. As the frequency is increased, a change in phase will be observed at the MMCs 
and the correlation changes. 

The phonon fluctuation noise between the MMC and the bath shows the opposite. At low frequency a 
positive correlation arises because the absorber has time to equilibrate and the phonon bath noise measured 
at MMC 1 will be sensed on MMC 2 in phase, so when combined the two contributions add. The positive, 
in-band, correlation of the phonon bath noise means that their contributions to the total summed noise are 
added linearly and not in quadrature. The phonon bath noise for a single pixel MMC with bath conductance 
G, will therefore be the same as a MMC DROID, with each MMC-bath conductance of G/2. 

At low frequency, the thermal fluctuations between the spin-electron system from MMC 2, measured on 
MMC 1 , tend towards zero. This is because, on long time scales the cold bath acts as a heat sink with infinite 
heat capacity; thus, in the zero frequency limit, noise fluctuations will flow to the heat bath and not be sensed 
by MMC 2. Consequently, at zero frequency there exists no correlation between the spin electron noise on 
each MMC. As the frequency is increased, more heat will flow through the absorber and a correlation arises. 


3.0 Energy and Position Resolution 


Under the regime of optimal filtering [17], the energy resolution of a microcalorimeter is calculated from the 
noise-equivalent-power ( NEP ) [18]. The energy resolution of a single pixel MMC, assuming just an electron 
and spin subsystem, can then be written as [14]: 


^Efwhm ~ 2.35^Jk b T C e 
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and is limited by fluctuations in energy exchange between the spins and the electrons. Here, To (~10‘ 6 s) is the 
relaxation time constant between the spins and the electrons and x\ (~10" 4 s) is the thermal time constant 
between the detector and the heat sink. The parameter /3 = C z /(C z +C e ), is the normalised heat capacity ratio. 
In the context of a canonical ensemble of thermal subsystems, as in our FEM, r 0 is defined by the thermal 
equilibrium time constant between the spin and electron system and thus depends explicitly on C z and G z . 
However, in reality the thermal conductance G z between the spin system and electron system is not an 
adjustable parameter. The time constant r 0 is more appropriately defined by the intrinsic, material dependent, 
spin-electron relaxation time (which varies inversely with temperature [14]) and as such is independent of 
the parameters C z and G z . Thus, for a single pixel MMC, with a fixed ratio t 0 /ti, Equation 8 is minimised 
when C z = C e (ft = _). Given the complexity of the DROID system, the appropriate optimisation of heat 
capacities is less obvious. We can however attempt to put some limits on what the optimisation should be. 
Assuming a MMC DROID with a short absorber thermal equilibration time r a bs (sufficiently large G a bs) 
compared to the thermal time constant between the MMC and heat sink T\, the MMCs will be strongly 
thermally coupled to the absorber electron system, and our detector will behave more like a single MMC. 
Then, from Equation 8, it follows that the detector energy resolution will be optimised when C z = C e + C a tJ2. 
However, in reality, G a bs must be designed to give enough spatial variation in pulse shape across the absorber 
to enable position determination. If r a b S and Ti are of similar time scales the thermal coupling of the MMCs to 
the absorber will be weaker and this optimisation may not be appropriate. In this regime, optimisation of the 
individual sensors, such that C z = C e , may be more applicable. 


Despite the limitations of our FEM in truly representing the physical nature of the MMC DROID, we can 
investigate the role thermal diffusion (determined by G a bs ) in the absorber has on energy and position 
resolution for a fixed set of MMC design parameters. The use of optimal filtering to numerically calculate 
energy and position resolution for position-sensitive TES DROIDs is described by ref. [12] and can be 
applied to the equivalent MMC concept, though we note these algorithms are yet to be implemented on real 
data. For this analysis we use the parameters in Table 1, which are estimates derived using the data in ref. 
[14], for optimised energy sensitivity of each individual sensor (assuming a typical single pixel Er doped Au 
MMC). 

For a TES based DROID, we have shown that the resolution depends explicitly on the ratio of the thermal 
conductances GabJG, and is independent of their exact values (though the permissible count rate is not) [12]. 
However, only in the physically unrealistic limit r 0 — » 0, will the MMC DROID depend only on G a bJG. 
Figure 5 shows the calculated energy resolution at the centre of the device, as a function of G a bs, evaluated 
for different G values, all other parameters are fixed. Only the intrinsic detector noise sources are included. 
Optimised energy resolution requires large G a bJG, that is, rapid thermal diffusion across the absorber. For a 
segmented absorber design (such as in [6]), the position resolution depends on the pixel pitch (typically 250 
pm is adequate for most astrophysics applications). However, for a continuous absorber design, the 
theoretical position resolution can be calculated from the optimal filtering algorithms described in ref. [12]. 
Figure 6 shows the position resolution at the centre of the device and Figure 7 shows the trade-off between 
energy and position resolution for the different G and G ahs . The position resolution is inversely proportional 
to incident photon energy and is calculated assuming a 1 keV photon. The minimum in position resolution 
occurs because of the requirement to have enough variation in the pulse shape across the absorber to enable 
position determination, while maintaining optimum noise performance. 

The error bars for the G = 1 nWK' 1 line in Figure 7, show the effect of including the SQUID flux noise and 
the magnetic Johnson noise. The effect on the energy resolution is less than 8 % for all G a bs- The effect on 


the position resolution is more significant, especially away from the minimum, where there is a ~ 20 % 
degradation in position resolution at G a b S = 100 nW K' 1 . These effects will be worse if the MMC sensitivity 
(dM/dT) is reduced, thus, particularly for high position resolution applications, optimisation of the device 
sensitivity is essential. 

A potential, disadvantage of the large absorber MMC is in the count rate. The decay time of the pulse r and 
hence the count rate, depends on the thermal conduction to the bath G. A MMC does not have the advantage 
of a reduced intrinsic (C/G) time constant as for a TES detector operated with electrothermal feedback 
(where typically r~ 100 ps) [1], Furthermore, the dependence of the time constants on resolution (Equation 
8) means there is an intrinsic trade-off between energy resolution and count rate. For G a bs> G, the pulse 
decay time for the MMC DROID is given by r = C lot /2G, where C tot is the total heat capacity including both 
MMCs and the absorber. A fast decay time of ~ 250 pis is achievable for G = 10 nW K' 1 , and although the 
theoretical energy and position resolution are not optimised they still give eV energy resolution and pm 
position resolution in the keV energy range, comparable to the TES DROID. The advantages of the MMC 
DROID over the TES counterpart are likely to be related to more practical aspects. The reproducibility of the 
TES transition temperature and shape between sensors is not yet at the level required for the detector 
responses to be perfectly symmetric [19]. Such variations between sensors complicate the operation of the 
device and ultimately can result in a loss of resolution. Furthermore, non-linearity in the detector response 
and TES saturation at high photon energy, particularly for events absorbed close to the TESs [6], can also 
result in a loss of resolution. However, for the MMCs, the intrinsic reproducible, linearity of dM/dT [15] 
means the detector response is much more linear with energy and the sensor responses should be symmetric. 
Furthermore, over the typical operating energy range the MMC will not suffer saturation at high photon 
energy as with a TES detector. These factors make the pursuit of a position-sensitive MMC an exciting 
prospect. 


4.0 Conclusions 


Using a finite-element model we have investigated the theoretic noise response of a position-sensitive, one- 
dimensional imaging spectrometer based around MMC technology. We have numerically calculated the 
theoretical best achievable energy and position resolution, assuming optimal filtering, across a 2 mm device 
using typical MMC parameters from literature. The analysis reveals the trade-offs in designing such 
detectors for different applications and provides a framework for the development of a practical device, 
which would further extend the possibilities for cryogenic focal plane instruments in X-ray astronomy. The 
practical advantages such as reproducibility of the signal response between sensors means that ultimately the 
development of MMC DROIDs may offer a technologically simpler alternative to the TES DROID. We have 
also extended the analysis to investigate two-dimensional systems such as a square absorber with MMC 
readout in each comer. In general, the noise correlations between sensors and the trade-offs between energy 
and position resolution follow the same scaling as for the one-dimensional equivalent. However, for a given 
number of readout sensors per absorber pixel, the total heat capacity will be twice as big compared to the 
one-dimensional case. Thus the energy resolution will be 2 m worse for an array of two dimensional 
DROIDs, compared an array of one dimensional DROIDs with the same number of readout channels. 
Consequently, the development of the two dimensional devices offers a less attractive alternative. 
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Table 1 . Parameters used to calculate the MMC DROID noise spectral density. 


Parameter 

Value 

c z 

IPJK' 1 

c e 

1 pJ K' 1 

Cabs 

lpJK" 

Cabs 

25 nW K’ 1 

G z 

500 nW K"' 

G 

1 nW K' 1 1 

L 

2 mm 

dM/dT 

44342 

T 

50 mK 

V 

2050 p,m 3 

r 

10.7 p,m 

Gk 

0.5 





Figure Captions 


Figure 1. Schematic diagram of the finite-element model used to calculate the response of the position-sensitive 
MMC concept. Only one MMC and part of the absorber is included. 

Figure 2. Flux noise spectral density (pd> 0 VHz) in the pickup coil of MMC 1 in a two MMC detector. Included 
is the phonon-bath noise (blue triangles) and spin-electron noise (red diamonds) associated with MMC 1, with 
spectra similar to that expected from a single pixel MMC [14], The phonon-bath (purple crosses) and spin- 
electron noise (orange pluses) from MMC 2 are all coupled through the absorber and therefore appear on the 
MMC 1 read-out. The thick dashed line gives the total noise contribution from each absorber element. The thick 
solid line gives the total noise spectrum. The thin solid line is signal response from a 1 keV photon absorbed at 
the centre of the device. 

Figure 3. Flux noise spectral density (pOo/VHz) when both MMC data is summed. The key is the same as in 
Figure 1. Due to symmetry of the device, some contributions are overlaid. 

Figure 4. Correlation coefficient, defined as the normalised covariance (<1>i<E> 2 (/))/( | <h n i (/) | )/( | ^n 2 (/)|>, 
between the noise measured on the two MMCs. Included is the total correlation (black solid), phonon-bath (blue 
dashed), spin-electron noise (red dot-dashed). 

Figure 5. Energy resolution vs. G a bs for four different values of G. 

Figure 6. Position resolution vs. G a bs for four different values of G. 

Figure 7. Trade-off between energy and position resolution for different G and G ahx . For the plot corresponding 
to G = 1 nW/K, the error bars indicate the resolution including the contributions from the SQUID noise and 


magnetic Johnson noise. 
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